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Abstract 


Tvo  procedures  for  specifying  element  weightings  for  antenna 
arrays  are  described.  The  first  procedure  is  used  for  gain 
optimization  and  the  second  for  shaping  the  array  gain  pattern. 
Applications  of  these  procedures  including  beamwidth  and  sidelobe 
reduction  are  included. 
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1.  INTRODUCTION 


This  note  presents  two  beamforming  algorithms  applicable  to 
arrays  of  antenna  elements.  Both  algorithms  are  derived  from  the 
standpoint  of  determining  a  set  of  system  parameters  to  optimize  a 
specified  system  performance  index  and  both  utilize  a  constrained 
optimization  approach. 

The  first  algorithm  addresses  the  problem  of  gain  optimiza¬ 
tion  for  arbitrary  arrays.  The  formalism  and  the  solution  of  this 
problem  is  an  illustration  of  the  minimum  energy  method  (1)  which 
is  closely  related  to  the  maximum  likelihood  method  (2»  3)  used  in 
adaptive  arrays.  It  is  presented  primarily  as  a  means  of  intro¬ 
ducing  a  formalism  which  can  be  applied  to  a  second  algorithm  that 
addresses  the  question  of  shaping  the  beam  to  achieve  a  desired 
antenna  pattern.  This  latter  algorithm  can  be  used  to  control 
both  the  amplitude  and  phase  of  the  antenna  pattern  at  specified 
angles  and  thereby  provides  a  solution  to  a  variety  of  problems 
faced  by  the  array  designer.  Examples  of  the  solutions  included 
in  this  paper  are  maximization  of  gain,  beamwidth  reduction  and 
control  of  sidelobe  levels.  Other  problems  which  are  amenable  to 
this  approach  include  not  only  amplitude  shaping  of  array  patterns 
but  also  phase  front  shaping. 

The  two  beamforming  algorithms  are  described  in  the  second 
section  of  this  note  followed  by  examples  in  the  third  section  and 
a  summary  in  the  final  section. 


2.  ALGORITHM  FORMALISM 
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Consider  a  linear  array  of  N  elements,  not  necessarily 
equally  spaced,  as  shown  in  Fig.  1.  Each  element  has  a  power 
pattern  v(6,  4)  a  complex  amplitude  of  excitation  a^,  1*1,..., N 
and  is  located  at  d^,  i~l,...,N.  The  array  factor  which  is  the 
electric  field  radiated  by  an  array  of  isotropic  elements  can  be 
written  as 

E(  0)  *  At«0*  (1) 

where  the  vectors  A  and  are  the  complex  column  vectors  defined 
below,  T  denotes  the  transpose  and  *  the  conjugate  operation. 


(2) 


Here,  k  is  the  wavenumber  defined  as  2lt/X,  where  X  is  the  wave¬ 
length  of  the  radiated  field.  The  total  power  density  radiated  by 
this  array  is  therefore  given  by: 

S(0,  ♦)  -  v(0,  *)  |  E(  0)  |  2 

where 
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Here  the  f  symbol  signifies  Hermitian  conjugate  and  l  is  the  rank 

0 

one  matrix  defined  below: 


^0  ^©t. 
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The  total  power  radiated  over  4w  steradians  is  given  by: 


p  "  4» 


which  can  be  written  as 


sin6S(8,  ♦)d8 


P  -  At  0A 


where  the  matrix  $  has  elements  $mn 


2*r  w- 

=  1  I  J  v(0,  ♦)  e“ik(dn“dm)cose  sine  de. 


In  order  to  maximize  the  gain  of  this  array  in  a  specific 
direction  (0O,  <fr0)  we  minimize  the  total  power  radiated  over  4ir 
steradians,  P,  subject  to  the  constraint  that  the  power  radiated 
in  the  direction  (60,  ♦  <,)  is  equal  to  a  constant,  for  example, 
unity.  We  therefore  minimize  the  total  power  radiated  over  all 
angles  except  in  the  desired  direction  (60,  ♦  „).  This  optimiza¬ 
tion  is  equivalent  to  maximizing  the  gain  S(60,  *oVp  in  that 
direction.  This  is  accomplished  by  using  the  technique  of 
Lagrange  multipliers  and  minimizing  the  cost  function 

J  *  P  +  X  [1-S(  e0,  ♦oH.  (7) 

For  an  isotropic  element  power  pattern,  (i.e.  v(e,  $)»1),  eq.  (7) 


reduces  to 


where 


J  »  At0A  +  X  [1-At£  A] 

6  o 

♦mn  “  sin  [k<dn“dm)l/[k(dn-dm)]. 


Taking  the  derivative  with  respect  to  A  and  equating  to 
zero  we  have 

3  J 

-  -  20A  -  2XT  A  -  0  (9) 
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resulting  in  the  eigenvalue  equation 

♦  A  =  X_1A  -  yA  (10) 

®  o 

The  eigenvalues  and  eigenvectors  of  0’1!  can  be  calculated  by  a 

®o 

simulataneous  diagonalization  of  l  and  0  but  the  solution  of 

®  o 

(10)  is  more  easily  obtained  by  recalling  that  £  is  a  matrix  of 

®o 

rank  one.  The  solution  for  the  optimum  weights  A,  is  the 

eigenvector  of  (10)  corresponding  to  the  single  non-zero 

eigenvalue.  It  is  straightforward  to  show  that  this  solution  is 

A  «  0-1  €0  (11) 

o 

and  the  non-zero  eigenvalue  is 

y  -  €te  0-1  €g  (12) 
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It  can  also  be  easily  shown  that  the  gain  of  the  array  at  60, 
i.e.  the  maximum  gain,  is  equal  to  the  eigenvalue,  y.  Equations 

(11)  and  (12)  represent  the  complete  solution  of  the  gain 
optimization  problem. 

We  next  turn  our  attention  to  an  algorithm  for  specifying 
the  complex  field  pattern  of  the  array  at  a  number  of  specified 
angles  while  minimizing  the  total  radiated  power.  Solutions  of 
this  problem  can  be  used  for  designing  arrays  with  specified 


beamwidths,  reducing  sidelobes  in  particular  directions,  altering 
phase  fronts,  etc.  The  formalism  for  this  problem  is  similar  to 
the  one  described  above,  but  in  this  case  the  constraints  are 
levied  on  the  radiated  fields  in  specified  directions  rather  than 
the  radiated  power  in  one  particular  direction.  The  constraints 
take  the  form 

E(0i)  =  ATe*  =  Ci  i=l, . . . ,ra  (13) 

®i 

where  c^,  the  amplitude  of  the  field  in  the  direction  0^  is  a 
complex  number  and  the  number  of  constraints,  m,  is  smaller  or 
equal  to  the  number  of  elements  in  the  array. 

The  criterion  function  to  be  minimized  is  given  by 

m 

J  *  At0A  +  l  Ui(Ci-AT£*  )  (14) 

i=l  ei 

where  y*  are  the  Lagrange  multipliers. 

The  solution  for  the  antenna  element  weights  is  given  by 

m 

A  -  1/2  l 

i=l  ei  (15) 

where  the  values  of  the  Lagrange  multipliers  y^  may  be  determined 
from  the  m  constraint  equations  shown  in  (13)  i.e. 

m 

l  Wj  0-1  €  -  2Ci  (16) 

j=l  J  6i 

Equations  (15)  and  (16)  completely  specify  the  desired  solu¬ 
tion  for  the  antenna  element  weights  for  the  beamshaping  problem. 
They  may  be  combined  into  a  single  equation  by  defining  the  matrix 


l.A  1..J 
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and  the  vectors  fi  andC 


so  that  (16)  may  be  written  as 


and  therefore 


hV^u  -  2C 


2(hV1h>-1C 


resulting  in  the  following  equation  for  the  antenna  element 


weights 


a  -  ♦_1h(hV1h)“1C  . 


(17) 


Equations  (11)  and  (17)  represent  the  major  result  of  this 
work.  The  former  prescribes  the  weights  of  the  antenna  elements 
for  maximizing  the  array  gain  in  a  particular  direction  and  the 
latter  prescribes  the  element  weights  for  shaping  the  beam  to  meet 
specific  constraints. 


3.  EXAMPLES 


In  this  section  we  apply  the  equations  of  the  previous 
section  to  the  problem  of  reducing  the  beamwidth  achievable  by  the 
conventional  endfire  and  broadside  arrays.  We  will  show  that  by 


application  of  eq.  (17)  the  beamwidth  can  be  reduced  at  the  ex¬ 
pense  of  a  reduction  in  gain  and  a  resulting  increase  in  sidelobe 
level.  For  the  case  of  the  endfire  array  we  include  an  example  of 
sidelobe  reduction  in  addition  to  beamwidth  reduction  to  further 
illustrate  the  application  of  eq.  (17). 

We  consider,  for  both  examples,  an  array  of  five  isotropic 
elements,  equally  spaced  at  5  inch  intervals  operating  at  2.0  GHz, 
which  results  in  a'0.85X  spacing  between  elements.  We  first 
examine  the  case  of  an  endfire  array,  that  is  the  element  weights 
and  phases  are  chosen  to  provide  maximum  gain  in  the  direction  of 
the  line  of  the  array. 

Figure  2a  shows  the  gain  pattern  of  the  array  obtained  from  a 
solution  of  eq.  (11)  for  maximizing  the  gain  in  the  direction 
9 o=0° .  A  similar  result  would  be  obtained  using  eq.  (17)  for  the 
case  of  a  single  constraint  in  the  direction  of  0i=O  deg.  This 
conventional  end-fire  array  solution  results  in  an  endfire  gain  of 
5.8  and  a  beamwidth  (measured  between  the  half-power  points)  of 
approximately  50  deg.  Using  eq.  (17)  to  reduce  the  beamwidth  to 
30°  we  obtain  the  gain  pattern  shown  in  fig.  2b.  As  a  result  of 
reducing  the  beamwidth  by  40%,  the  endfire  gain  has  been  reduced 
by  45%  to  a  value  of  3.15  and  the  total  power  radiated  increased 
by  80%.  In  addition,  the  gain  in  the  broadside  direction  has 
become  appreciably  greater  than  in  the  endfire  direction.  In  many 
applications,  this  particular  sidelobe  level  increment  may  not  be 


particularly  worrisome  but  there  may  be  some  applications  where  a 
large  gain  at  90  deg.  may  bo  undesirable.  Figure  2c  shows  the 
result  of  reducing  the  gain  at  90  deg.  while  maintaining  the 
endfire  beamwidth  at  30  deg.  The  endfire  gain  has  now  been 
reduced  to  a  value  of  2.4,  which  is  less  than  half  of  the  original 
gain  shown  in  fig.  2a  and  the  total  radiated  power  is  2.4  times 
greater  than  the  value  shown  in  that  case.  We  have,  however,  been 
able  to  achieve  a  30  deg.  beamwidth  and  in  addition,  significantly 
reduce  the  extremely  high  broadside  sidelobe  level  incurred  in  the 
intermediate  result  in  which  beamwidth  control  was  the  only 
constraint.  The  cost  of  this  specific  sidelobe  reduction  has  been 
an  increase  in  overall  sidelobe  level  over  the  entire  angular 
range. 

Turning  to  the  case  of  a  broadside  array,  we  see  the  result 
of  maximizing  the  gain  in  fig.  3a.  The  beamwidth  of  the  broadside 
lobe  is  15  deg.  for  this  array.  Reducing  the  beamwidth  to  9° 
results  in  the  gain  pattern  indicated  in  fig.  3b.  The  broadside 
gain  has  been  reduced  by  45%  and  the  radiated  power  increased  by 
80%  as  a  result  of  this  40%  reduction  in  beamwidth  . 
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(a)  GAIN  OPTIMIZATION 


Fig.  2.  Endfire  Array  Gain  Patterns 


This  result  is  similar  to  that  achieved  for  the  case  of  the 
endfire  array  in  that  the  reduction  in  mainlobe  gain  and  increase 
in  radiated  power  level  are  the  same  for  an  equivalent  amount  of 
beamwidth  reduction.  The  increase  in  sidelobe  level  is  evident  in 
the  increased  gains  at  70  and  110  deg.  as  well  as  at  bther  regions 
of  the  gain  pattern.  However,  in  this  case  unlike  the  endfire 
case  there  is  no  large  increase  in  a  particular  sidelobe  so  we 
will  not  illustrate  sidelobe  reduction  for  this  array. 

4.  SUMMARY 

A  procedure  for  determining  the  weights  of  array  elements  to 
shape  both  the  amplitude  and  phase  of  the  radiated  field  has  been 
presented  along  with  an  algorithm  for  determining  the  weights  that 
maximize  the  gain  in  a  specific  direction. 

The  beamshaping  procedure  is  a  constrained  optimization 
method  which  minimizes  the  total  radiated  power  while  constraining 
the  complex  value  of  the  radiated  field  at  specifc  angles.  Vhile 
this  provides  both  amplitude  and  phase  control  over  certain 
portions  of  the  radiated  field,  the  value  of  the  gain  cannot  be 
controlled  over  these  regions.  This  results  in  trade-offs  between 
array  parameters  (e.g.  gain,  beamwidth  and  sidelobe  level)  that 
the  array  designer  must  judge  for  each  application.  Examples  of 
typical  results  and  trade-offs  were  presented  for  the  case  of  the 
conventional  end-fire  and  broadside  arrays. 


REFERENCES 


1.  D.  H.  Johnson  and  S.  R.  DeGraaf,  "Improving  the  Resolution  of 
Bearing  in  Passive  Sonar  Arrays  by  Eigenvalue  Analysis,"  IEEE 
Trans.  Acoust.,  Signal  Processing,  Vol.  ASSP-30,  pp.  638-647, 
1982. 

2.  J.  Capon ,  "Maximum-Likelihood  Spectral  Estimation,"  Nonlinear 
Methods  of  Spectral  Analysis,  Topics  in  Applied  Physics, 


Volume  34,  ed.  by  S.  Haykin  (Springer-Verlag,  Nev  York,  1979). 

3.  R.  N.  McDonough,  "Application  of  the  Maximum-Likelihood  Method 
and  the  Maximum- Entropy  Method  to  Array  Processing,"  Nonlinear 
Methods  of  Spectral  Analysis,  Topics  in  Applied  Physics, 

Volume  34 t  ed.  by  S.  Haykin  (Springer-Verlag,  New  York,  1979). 


13 


MCUMTY  ClAMtnCATKMi  Of  THIS  PAGE  (Whom  Bum  bM| 


32 MESm 


ESD-TR-82-156 


4.  TTTli  (WW  Subtitle) 


A  Note  on  Beamforming  Algorithm* 


7.AUTH0MV 


William  H.  Schoendorf 


Technical  Report 


7 . I  •  l  rrr  ,i 


Technical  Report  631 


i  cogtract  on  mast  mmew 

F 19628- 80- C-0002 


Lincoln  Laboratory,  M.I.T. 
P.O.  Box  73 

Lexington,  MA  02173-0073 


OFFICE  MMEASOAMMESS 

Air  Force  Systems  Command,  USAF 
Andrews  AFB 
Washington,  DC  20331 


M  MOMTOMM  ASEICY  MME  a  A00MSS  (if  different  from  Controlling  Office) 


Electronic  Systems  Division 
Hanscom  AFB,  MA  01731 


Project  No.  627A 
Program  Element  No.  6331 IF 


tl  REPOST  DATE 

12  January  1983 


IL  MCUMTY  CIASS.  (of  this  report) 
Unclassified 


IS.  MSTMBtmOS  ST4THKST  (of  tkU  Report) 


Approved  for  public  release;  distribution  unlimited. 


17.  MfTMStfTNM  STATEMEST  (of  the  abetraa  emend  in  Bloch  BO,  if  different  from  Report) 


rnrTCTrn  ttitth 


II.  KEY  W0MI  (Continue  on  reveree  aide  if  necemnry  and  identify  by  Hock  number) 


antenna  arrays 


gain  optimization 


array  gain  pattern 


IS.  AM  TRACT  ( Continue  on  reveree  tide  if  neremary  on d  Identify  by  block  number) 


Two  procedures  for  specifying  element  weightings  for  antenna  arrays  are  described.  The  first  pro¬ 
cedure  is  used  for  gain  optimization  and  the  second  for  shaping  the  array  gain  pattern.  Applications 
of  these  procedures  including  beamwidth  and  zidelobe  reduction  are  included.  . 


00  FMM  1473  EDtTKM  OF  1  MV  M  IS  08SOUTI 


UNCLASSIFIED 


I  Jm7> 


MCUMTY  CLASSIFICATION  OF  THIS  PAGE  (When  Deem  Soured) 


